Electrostatic interactions of charged bodies from the weak to the strong coupling regime 
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A simple field theory approach is developed to model the properties of charged, dielectric bodies and their 
associated counterions. This predictive theory is able to accurately describe the properties of systems (as com- 
pared to computer simulation data) from the weak coupling limit, where the Poisson-Boltzmann theory works 
well, through to the strong coupling limit. In particular, it is able to quantitatively describe the attraction between 
like-charged plates and the influence of image charge interactions. 



Electrostatic interactions play a major role in determining 
the structure and thermodynamics of many colloidal and bi- 
ological solutions, which typically contain charged macro- 
molecular structures with low dielectric interiors, such as 
DNA, charged micelles, or membranes. These charged struc- 
tures are always surrounded by neutralizing counterions, and, 
in many cases, the properties of the system can be mainly at- 
tributed to properties of the counterions Jl]] . 

When the electrostatic interactions are weak, their contri- 
butions to the system properties are accurately described by 
the Poisson-Boltzmann (PB) theory. However, as these inter- 
actions strengthen, the PB theory becomes less and less accu- 
rate. Perturbation methods, such as the loop expansion can be 
used to systematically improve the theory; however, the first 
loop correction offers only a small improvement JUIUS], and 
higher-order corrections are increasingly complicated to eval- 
uate. 

When the electrostatic interactions are extremely strong 
(e.g., when the surface charge of the macromolecular struc- 
tures or the valency of the counterions is high), the Poisson- 
Boltzmann theory can yield qualitatively incorrect predic- 
tions. For example, in this regime, the counterions can gen- 
erate attractive forces between similarly charged objects ]4[]. 
This phenomenon cannot be explained by the PB theory, but 
has been observed in Monte Carlo simulations J4[] and in ex- 
periments (e.g., condensation of DNA molecules [5], bundle 
formation of filamentous actin |6]. In this strong coupling 
regime, the counterions "collapse" on the neutralizing charged 
surfaces to form a highly interacting 2D structure that resem- 
bles a confined one-component plasma (OCP) J3, B lal- In 
these systems, the average distance a± between the ions is 
much larger than the average distance z between the ions and 
the charged surface (i.e. z -C a±). Consequently, a single par- 
ticle theory provides a good description of the system. This 
leads to the strong coupl ing (SC ) expansion Jf! [Toll , which 
has been quite successful ll llll2lll3l1 . 

However, many systems are in a regime where both the PB 
theory and SC expansion are inaccurate. The behavior of these 
systems can be rationalized in terms of a correlation hole Um 
— a region of size a around each counterion where it is unfa- 



vorable for other counterions to be located. At length scales 
greater than a, the counterions are weakly correlated, while 
at shorter length scales, the counterions are strongly corre- 
lated but fairly "isolated" 1 1011 . In the weak coupling regime, 
the counterions form a diffuse 3D layer, and the size of the 
correlation hole is approximately equal to the Bjerrum length 
Ib = ftq/e (where q is the counterion charge), the distance 
at which two counterions interact with energy ksT. In the 
SC limit, the size of the correlation hole becomes equal to the 
average (2D) distance between the ions aj_ = 2il^/21b Jp>, 
where /i = (2ir(3q'E<)~ 1 is the Gouy-Chapman length (where 
S is the surface charge density), the distance at which the in- 
teraction between a counterion and the charged surface equals 
k B T. 

Based on this observation, Weeks and coworkers lfl5ll and 
Santangelo [16] developed approaches that split the interac- 
tion between the ions at short and long range. The long- 
range interaction is treated within a mean field approximation, 
and the short-range interactions with a more precise approach 
(e.g., computer simulation, liquid state theory, etc.). With an 
appropriate value for a, these approaches can successfully de- 
scribe Monte Carlo results for the full range of electrostatic 
coupling. However, the value of a is determined empirically. 
Additionally, these approaches are not capable of describing 
systems with dielectric inhomogeneities. 

In this work, we present a self-consistent theory that is in 
good agreement with Monte Carlo simulations at weak, inter- 
mediate, and strong coupling. The theory is similar in idea to 
the work of Weeks and coworkers [15] and Santangelo lllrjll . 
however, the parameter a is calculated consistently from the 
partition function, rather than chosen empirically or adjusted 
to fit data. In addition, the theory accurately describes the 
presence of dielectric bodies, even in the SC limit, which has 
not been demonstrated by any previous theory. For the two 
plate system in the presence of image charges, the system un- 
dergoes a transition from from a two peak density profile to a 
one peak density profile. 

We limit our attention to systems composed of a fixed 
charge distribution E(r) that is surrounded by a neutralizing 
cloud of counterions, which are point charges of magnitude 
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q immersed in a medium with dielectric constant e and pos- 
sibly in the presence of dielectric inhomogeneities. The total 
electrostatic energy H of the system is given by 



H=\j drdr'Q(r)G (r,r')Q(r') 

q 2 

- -^-Gf rco (r fc ,r fc ) 



(1) 



where Q(r) = (5 d (r — rfc) + S(r) is the total charge den- 
sity, rfc is the position of counterion k, Go is the Green's func- 
tion of the associated electrostatics problem (including the ef- 
fects of dielectric objects), and Gf ree is the Green's function 
in the absence of dielectric inhomogeneities. 

To se para te short and long wavelength phenomena, we split 
1 15, 17 1 the Green's function Go into a short wavelength 
G s and a long wavelength Gi component 



Go(r,r')=G s (r,r') + G;(r,r') 



(2) 



where G s = (1 - V)G , and G; = VG . The operator V 
filters out the short wavelengths; its specific form is arbitrary, 
and in this work we use V = [1 — <t 2 V 2 + (T 4 V 4 ] _1 , where 
a is the length scale which divides the long from the short 
wavelength phenomena. 

Equation (Q]i can be written as: 



H = \j dvdr'Q{v)Gi{Y,Y')Q{v') + E s 



E 



o 2 q 2 
u ( r k) ~ yG s (r fc ,r fc ) - —V5G {r k ,T k ) 



(3) 



where SGq = Go — Gf ree , u(r) is a one-particle interaction 
potential given by 



u(r) = q J dr'G s (r, r')S(r') + ^SG (r, r) 



(4) 



-VG{ iee (r,r), 



and E sc is the self energy of the fixed charges, defined as 

E SC = \J drdr'E(r)G s (r, r')E(r'). (5) 

By performing a Hubbard-Stratonovich transformation 
1911 twice on the grand partition function of the system, two 
fields are introduced: ipi, which is associated with G; and 
represents interactions at length scales greater than a, and 
tp s , which is associated with G s and represent interactions at 
length scales less than a. In the approximation scheme we 
pursue, the one-particle contribution to the partition function 
is treated exactly, while the interaction between the particles is 



treated approximately. The field ijj s is strongly fluctuating and 
coupled not only to itself, but also to the field tpi . To evaluate 
the functional integration over tp s , we use a cumulant expan- 
sion, truncated at first order. This leads to an approximation 
similar to the SC expansion of Moreira and Netz BlOluOll . 

Performing the functional integration over ip s , we get an 
effective field theory for the long wavelength system. The 
expression for the grand partition function Zq is 



InZ G [7,E] = -l3E sc + ln 



(6) 



and Hi is the effective Hamiltonian, which is a functional of 
the field tpi given by 

-H^i] « J drdv'i, l (r)G^ 1 (r,r'W l (r') 

+ [ rf r A- 3 e -T-« i ^«-' 3u ( r )+T^' 5G '( r ' r ) (7) 

drT,(r)itpi (r). 



The long-wavelength field ipi is weakly fluctuating, so a 
mean-field approximation is sufficient to evaluate the func- 
tional integral over tpi in Eq. ©: 



lnZ G [ 7 ,E] / drdr'i^(r)Gr 1 (r,r'l/V / (ri 



drS(r)^(r) - f3E sc 



(8) 



The value of the mean field ipi (r) is determined by solving 
the Poisson equation 
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V • e(r)V0(r) = S(r) + ^(r) 



(9) 



where the electric potential is defined as /3<p = V 1 iipt, and 
the counterion density distribution p(r) is given by 



(10) 



Note that the density depends on the field ■0/ =V(j> rather than 
on the electric potential <fi, as in the PB theory. In the systems 
we consider, where there are only counterions, the chemical 
potential 7 is determined by the electroneutrality constraint. 

All properties of the system should be independent of 
the parameter cr; however, because the theory is approxi- 
mate, there will be a dependence on the value of a. Based 
on this, we determine the value of a by requiring that the 
grand partition function is stationary with respect to cr (i.e. 
d In Zq [7]/ da = 0). This is similar to the optimized random 
phase approximation 112 ill . 

Now we apply the theory to a system of counterions con- 
fined to one side of a plate with dielectric constant e' and a 
uniform surface charge S(r) = <5(z)£, where z is the distance 
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FIG. 1: (a) Counterion density profile near a single charged plate 
with A = 0. Counterion density profile for (b) 3 = 10 and (c) 
5 = 1000 near a single charged, dielectric plate, (d) Dependence 
of the parameter a with the cou pling p arameter H. The symbols 
are Monte Carlo simulation data 1 11, 22, 23j], the thin lines are the 
predictions of the SC expansion fill, 123 , 12311 , and the thick lines are 
from the present work. 



from the surface of the plate. For this system, the one-body 
potential of the counterions (see Eq. reduces to 



0u(v) = -^-(l + A)e 



✓3; 



COS ■ 
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(11) 



where z is the distance from the surface of the plate, and 
A = (e' — e)/(e' + e). The self energy of the surface charge 
(see Eq. ©) is (3E SC /N = ^-(1 + A) where N is the 
total number of counterions in the system. The strength of 
the electrostatic interactions is characterized by the parameter 
5 = Ib/i-I- 

The deviations of the counterion density profile from the 
predictions of the PB theory for the single plate system are 
plotted in Fig. [T^ in the case when there is no dielectric in- 
terface (i.e. A = 0). As the strength of the electrostatic in- 
teractions increases, the PB theory overpredicts the repulsion 
between the counterions and the ion cloud associated with the 
charged surface, and, consequently, it underestimates the ad- 
sorption of the counterions to the surface. The SC limit (given 
by the thin line) is only approached at fairly high values of 
the coupling parameter (H > 100). The predictions of the 
present theory are, however, in good agreement with Monte 
Carlo simulation data, exhibiting the crossover of the density 
profile from the PB theory to the SC limit with increasing val- 
ues of 2. 

When the plate has a low dielectric interior (i.e. e' < e, 
A > 0), repulsive image charge interactions with the plate re- 
pel the counterions from the surface, which oppose the attrac- 
tive interactions with the surface charge. Consequently, the 



counterions no longer collapse onto the surface, but instead 
peak at a distance away from the surface |24, 25, 2^1 . Coun- 
terion density profiles for this situation are plotted in Figs.[TJ) 
and c. When A > 0, the SC expansion fails to accurately 
describe the Monte Carlo data, even well into strong coupling 
regime (H = 1000). The failure of the SC theory is due to 
the fact that the average distance between the ions and the dis- 
tance to the charged plate are the same order of magnitude. 
Interactions between the counterions become significant, and 
the physical basis for the SC is no longer fulfilled. 

Our approach is able to overcome these difficulties because 
in addition to including the one-body interactions with the 
plate interactions, which occur at length scales less than a, it 
also accounts for the interaction between ions at length scales 
greater than a. The variation of a with the coupling parameter 
S is plotted for various values of A in Fig. [TJl. For high val- 
ues of the coupling parameter, the size of the correlation hole 
quickly approaches a constant with respect to the spacing a± 
between the counterions. As A increases, the size of the cor- 
relation hole decreases, reflecting the increasing importance 
of the interactions between the counterions. 

Now we consider systems where the counterions are con- 
fined between two uniformly charged plates, separated by a 
distance d. Fig. [2^ shows the counterion density profile where 
the plates have dielectric constant e' = e and are separated 
by a distance d//i = 2; results for the pressure are shown 
in Fig. I2J3. At weak couplings, the force between the plates 
is strictly repulsive, but as 3 increases, a region of attraction 
develops at intermediate plate separations. The predictions 
of this work are in quantitative agreement with the computer 
simulation data 1 10], while the SC limit is applicable only at 
extremely high values of the coupling parameters (5 ~ 10 5 
for the pressure, see Ref. ifioll ). 

Decreasing the dielectric constant of the plates (e' < e) 
leads |fl2[|27ll to a qualitatively different counterion density 
profile, as shown in Fig. [3^. For sufficiently large plate sep- 
arations, there is a counterion peak next to each plate. When 
the distance between the plates is smaller or comparable to 
the average distance between the ions (aj_//i w \/2S, so 
d/fi < V2S), the repulsive image charge interactions push 
the counterions into a single peak in the middle of the plates. 
At these separations, the distance between an ion and its im- 
age charge is comparable to the average distance between the 
ions. The variation of the mid-point counterion density with 
the coupling constant is shown in Fig.fSJ?. The predictions of 
the present work are in good agreement recent Monte Carlo 
simulation results for all conditions examined. The SC theory 
II2II13J] is only accurate when d/fi <C V23. The parameter 
a (see Fig. [3}:) increases with increasing values of EE and de- 
creasing values of A, in a similar manner with the one plate 
case. 

The approach developed here is applicable to general ge- 
ometries (e.g., spherical or cylindrical) in the same manner as 
the Poisson-Boltzmann equation. The only difference will be 
the form of the fixed charge density S(r). Once this is given, 
the one particle potential is obtained by Eq. and the den- 
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FIG. 2: (a) Density profile for counterions confined between two 
charged plates with A = separated by distance d = 2/i. (b) Force 
between the plates as a function of their separation d. The symbols 
are simulations data [ 10], the thin lines are the prediction of the SC 
theory 1 10], and the thick lines are from this work. 



sity by Eq. (fTOt . The electric potential is obtained from the 
solution of the Poisson equation (see Eq. (0), which can be 
solved by standard methods fislllql . This theory can also be 
systematically improved by either increasing the order of the 
cumulant expansion for the integration over tp s or by going 
beyond the mean field theory for ipi . 

We thank Prof. C. Holm and Dr. M. Sega for providing us 
with unpublished simulation data for the two plate system. 
MM Hatlo acknowledges support from an EC Marie Curie 
Fellowship (MEST-CT-2004-503750). 




FIG. 3: (a) Density profile for counterions confined between charged 
plates with A = 0.98 and d = 10[i. (b) Counterion density at the 
mid-point between charged plates separated by distance d — 14[i 
and (c) the corresponding variation of the parameter a. The symbols 
are Monte Carlo simulation data 11311 . the thin lines are the prediction 
of the SC theory 1 10, 12, 13], and the thick lines are from this work. 
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